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Abstract 

The inclusive production of the C(;(782) vector meson in hadronic Z decays is mea- 
sured and compared to model predictions. The analysis is based on 4 million 
hadronic Z decays recorded by the ALEPH detector between 1991 and 1995. The 
production rate for Xp = Pmeson/Pbeam > 0.05 is measured in the uj —>■ 7r"*"7r~7r'^ de- 
cay mode and found to be 0.585 it 0.019stat i 0.033sys per event. Inclusive r] meson 
production is also measured in the same decay channel for Xp > 0.10, obtaining 
0.355 lb O.Ollstat i 0.024sys per event. The branching ratio for uj ^^/^^ is investi- 
gated. A total of 18.1 ±5.9 events are observed, from which the muonic branching ra- 
tio is measured for the first time to be BR(a; — > z^"*"/^") = (9.0±2.9statil-lsys) x 10~^. 
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1 Introduction 



The description of the hadronization process in QCD is deeply connected with the 
confinement property and requires non-perturbative methods. The precise measurement 
of momentum spectra of identified particles in the clean environment of e"'"e~ annihilation 
into hadrons may improve the understanding of hadronization. Meanwhile these 
measurements are necessary to test and tune the phenomenological models used to 
describe the hadronization: each of these models has free parameters which must be 
determined from comparison with data 

Resonant states and their dynamics are more closely related to the original partons 
than pions. The t] and ci;(782) mesons are copiously produced in hadronic events and thus 
well suited to a study of hadronization. 

In this investigation a measurement of the inclusive momentum distributions of the rj 
and u mesons obtained through their tt+tt^tt^ decays is presented, which improves on the 
previous ALEPH results p[ with higher statistics and reduced systematic effects. It also 
complements the recent publication concerning the 77 decay of the r] meson Studies 
by the OPAL and L3 Collaborations can be found in ||, 

The measurement of the partial width of leptonic decay rates of the light vector mesons 
is a good test both for the quark assignments (J2 a^qq with q = u,d, s) in the vector meson 
and the quark charges. With p° = {nn — dd) / V2, to = {nu + dd) / ^/2 and = ss the 
partial decay width of a vector meson V to lepton pairs (£=e^, /i^) can be calculated 
using the Van Royen-Weisskopf formula 0: 

where = (J2 o-qQqY is the square of the flavour-weighted sum of the charges Qq of the 
quarks (g = u, d, s), 'ipifi) is the amplitude of the qq wavefunction at the origin, and my is 
the meson mass. The expected ratio of the leptonic widths is Ti{p^) : Ti{uj) : T£{(f)) = 9 : 
1 : 2 in agreement with the widths and branching ratios measured as summarized in |^ 
for the decay to e^e~. While Ti{p^) and T£{(f)) are in agreement with the prediction based 
on the Van Royen-Weisskopf formula and lepton universality also for the muonic decay, 
only an upper limit so far has been established for the decay u (Table |l]). 

The muonic branching ratio of the u is measured here for the flrst time. 



Table 1: Widths and leptonic branching ratios of the light vector mesons V 





r [McV] 


BR(V-^e+e-) 


BR(V ^ /i+/i-) 




150.2±0.8 


(4.49 ±0.22) X 10-5 


(4.60 ±0.28) X 10-5 


UJ 


8.44±0.09 


(7.07 ±0.19) X 10-5 


(CL=90%) < 1.8 X 10-^ 




4.46±0.03 


(2.91 ±0.07) X 10"^ 


(3.7 ±0.5) X 10"^ 



After a brief description of the ALEPH detector, the selection of hadronic events 
is detailed in Section 3 and the particle identiflcation in Section 4. The improved 
measurement of inclusive r] and uj production in the tt+tt-tt'' channel with a large data 
sample of about 4 million hadronic Z decays is presented in Section |5|. The measurement 
of the branching ratio uo fi~^fi~ is presented in Section |^, followed by a summary and 
conclusions. 



2 The ALEPH Detector 



The ALEPH detector is described in detail elsewhere [y, Charged particles are 

measured over the polar angle range |cos^| < 0.85 and |cos^^| < 0.69 by the two layers 
of the silicon vertex detector (VDET). This is surrounded by a cylindrical inner drift 
chamber, and a large cylindrical time projection chamber (TPC) which measures up to 

21 three-dimensional space points per track. A particle's energy loss is sampled in the TPC 
by up to 338 wires and 21 pads. The tracking detectors are immersed in a magnetic field 
of 1.5 T and provide a momentum resolution of 5pt/vt = 0.0006pt © 0.005 {pt in GeV/c). 
The TPC is surrounded by an electromagnetic calorimeter (ECAL) of lead-proportional 
tube construction, which covers the angular range |cos^| < 0.98 and has a thickness of 

22 radiation lengths. It is finely segmented in projective towers of approximately 0.9° by 
0.9° providing an angular resolution of uq^^ = 2.5/v^ -|- 0.25 {E in GeV; ae,4, in mrad). 
The energy resolution is (Je/E = 0.18/v^ + 0.009 {E in GeV) for isolated showers. 
The hadron calorimeter (HCAL) uses the iron return yoke as an absorber, for a total of 
7.5 hadronic interaction lengths. The iron is interleaved by 23 layers of streamer tubes, 
which provide a two-dimensional measurement of muon tracks and a view of the hadronic 
shower development. The HCAL is used in conjunction with the muon chambers, which 
are two double-layers of streamer tubes with three-dimensional readout, and the tracking 
detectors to identify muons. The calorimeters and the muon chambers cover nearly the 
entire 47r solid angle. 



3 Event Selection 

For the event selection, good tracks are defined as originating close to the interaction 
point (with transverse impact parameter |(io| < 2 cm and longitudinal impact parameter 
l^ol < 5 cm), having at least 4 TPC hits, a polar angle in the range 20° < 9 < 160°, and 
a transverse momentum pt > 200MeV/c. Four million hadronic Z decays are selected by 
requiring at least five good tracks. The total energy carried by all good tracks is required 
to exceed 15 GeV and the sphericity axis must be in the range 35° < 6 < 145°. With 
these cuts, a sample of 3.0 million events is selected. The background to these events 



arises from tau pairs and two-photon events and is estimated to be less than 0.4% |jTT . 

For the measurement of the branching ratio BR(ci; — > /i^/i~) additional event cuts 
are applied to reject heavy flavour events and events with missing energy: events from 
Z — > bb and Z — >• cc are rejected by a lifetime tag [l^ keeping about 95% of light flavours. 



A further reduction is obtained by a cut on the missing energy corrected with hemisphere 
masses 



131 , -Emiss < 6 GeV, in the u hemisphere. 



For the purpose of comparing with models and as a means of measuring the detector 



acceptance, samples of events generated with the Jetset 7.4 Monte Carlo [|T4|, modified 
with DymuS [0 for electromagnetic radiative effects and improved bottom and charm 
decay tables, were passed through a full detector simulation and reconstruction program. 
The generator was tuned to describe the ALEPH data using the inclusive charged particle 



and event shape distributions [0. For model comparison and as a systematic check 
for the extrapolation into the unmeasured region the measured spectra were compared 
to those of Herwig 5.9 WTj. The model parameters were tuned using ALEPH data 



in the same manner as mentioned above. For the measurement of the branching ratio 
BR(co' f^~^fJ'~) a sample of events was generated with the Jetset 7.4 Monte Carlo, but 
with the muonic branching ratio of the u meson BR(co' — /U+/i~) set to 100%. After a full 



detector simulation and reconstruction, the generated events were subjected to the same 
selection and analysis chain as the data. 



4 Particle Identification 

The 7] and u mesons are produced and decay at or near the primary interaction point. To 
reduce background from charged particles produced far from the interaction point, e.g., 
from Kg decays or photon conversions, track candidates are required to be reconstructed as 
good tracks (as discussed above) with tighter cuts on the impact parameters {\do\ < 0.5 cm 
and l^ol < 3.0 cm) and the transverse momentum {pt > 250MeV/c). 

For the measurement of the branching ratio BR(u; —>■ fi^fJ^^) additional cuts are 
applied to identify and select muons above the background of mainly charged pions. 
One VDET hit is required. For the identification of muons the measured energy 
loss in the TPC, dE/dx^eas should be consistent with the expectation for a muon, 
(di?/dxmeas — d-E/dXgxp ) /c"exp > ~2, where dE/dXexp and cXexp are the expected energy 
loss and its resolution. The pattern in HCAL and the muon chambers should match with 
the pattern expected for penetrating muons, if the momentum is larger than 2.4GeV/c. 
For lower momenta, 1.5GeV/c < p < 2.4GeV/c, the number of fired planes within the 
last ten HCAL planes divided by the number of all fired planes should be between 0.35 
and 0.75 



For the selection of neutral pions two-photon invariant mass spectra are formed. The 
photon energy is estimated from the energy collected in the four central ECAL towers of a 
cluster, correcting to the full energy from a parametrization of the shower shape for a single 
photon in the calorimeter. While the energy resolution is degraded to a^i/i? ~ 0.25/a/E' 
{E in GeV) with this technique, hadronic background and clustering effects are reduced. 
Photon candidates are only accepted if the estimated energy is greater than 0.8 GeV. 
Neutral pion candidates are accepted if the energy of the photon pair is less than 16 GeV 
(to remove large uncertainties in the acceptance correction at high energy). The invariant 
mass of the photon pair must also be within ±3cr of the expected mass. The vr" energy 
resolution is improved by constraining the mass of the 7r° candidates to 135 MeV/c^. 
The poor purity at low momentum, due to the large multiplicity of low energy photons 
giving rise to a large combinatorial background, is improved by a "ranking" method: 
candidates that share photons with other candidates are ranked in an order determined 
by a estimator. The estimator R is calculated for each candidate from the photon pair 
opening angle 612 and the from the mass constraint: 

i? = ^12(1.0 + 0.1x'). (1) 

All 7r° candidates sharing photons with other candidates are removed except the one with 
the smallest value of R. 



5 Inclusive Production of r/ and u; 

5.1 Extraction of the 77 and a; Signals 

The rj and uj cross sections are extracted from invariant mass distributions in the tt+tt^tt^ 
decay mode. All charged tracks are assigned the charged pion mass; for neutral pion 
candidates the kinematic refit constrains the mass to the vr" mass. Data are analysed 
in six intervals of scaled momentum Xp = Pmeson/Pbeam, where Pbeam is the LEP beam 
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Figure 1 : Fits in the r/ region made to the invariant mass spectra of data (like-sign subtracted) . 

momentum. Only candidates satisfying Xp > 0.10 for the 1] and Xp > 0.05 for the a; are 
considered, as the signal to background ratio is too small for lower momenta. 

The invariant mass spectra are fitted as the sum of a background and a signal 
function. Subtracting the like-sign spectra from the unlike-sign spectra allows the 
backgrounds to be represented simply by a third order polynomial. The signal function 
is empirically parametrized as the sum of three Gaussian functions; the relative widths 
and normalizations of the three Gaussians are determined from the Monte Carlo and are 
fixed with the overall normalization left free. The results of fits to the data for the r] and 
uj are shown in Figs. |] and 0, respectively. The fit to the first momentum interval for the 
7] yields a large (~50%) uncertainty, therefore this momentum interval is not considered 
further. 



5.2 Measured Rates and Differential Cross-Section 

The production rate per event, R, is calculated for each momentum interval by correcting 
the fitted signal S for the reconstruction efficiency e and normalizing to one event. 
The rate is corrected for the branching ratio BR(?7 — > tt+tt^tt") = 0.230 ± 0.004, 
BR(cu ^ tt+vt-ttO) = 0.888 ± 0.007, and BR(7r° ^ 77) = 0.9880 ± 0.0003 @] to give 
the rate for all decay modes. The calculation is as follows: 

where N is the total number of hadronic events in data and e is the efficiency. is 
obtained as = A^gen ■ ^acc/^acc"; where Ngen is the number of generated events in the 
Monte Carlo, while N^^^ and A*"^^ are the numbers of accepted events after preselection 
in data and Monte Carlo, respectively. The efficiency e is defined as e = "^^rcc/'^'gen, where 
rircc is the number of reconstructed rj (a;) mesons in the Monte Carlo, and ngen is the 
number of generated rj [uj) mesons in the Monte Carlo (before event selection cuts). 
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Figure 2: Fits in the to region made to the invariant mass spectra of data (hke-sign subtracted). 



The implementation of the rj tt+tt^tt'^ decay in Jetset is according to phase space, 
while for the u tt+tt^tt'^ decay it is according to the correct matrix element. For this 
analysis the Monte Carlo was reweighted for the t] decay by parametrizing the decay 
transition probability as oc 1 — Tq/Tq^^^ as suggested in where Tq* is the kinetic 
energy of the neutral pion in the t] rest frame and T^^^^ is its maximum possible value. 
This form is consistent with A,, oc 1 + 2ay proposed in |]20[ and the measurements in 
2^, |2l|], where y = {3Tq/Q) — 1, Q is the mass difference between the r] and its daughters, 
and a = —0.47 ± 0.04. The t] efficiency is taken as the average of the results obtained 
with the two parametrizations. 

The overall efficiency in the measured momentum region is 15.3% for the rj and 11.2% 
for the uj. The measurements cover only the momentum region Xp > 0.10 for the rj 
and Xp > 0.05 for the uj. To estimate the total production rate the measured rates are 
extrapolated to a;p = using the fragmentation function in the Monte Carlo. The fraction 
of rj mesons for Xp > 0.10 is 29.56%. For the uj meson the fraction for Xp > 0.05 is 58.8%. 



5.3 Systematic Uncertainties 

The largest systematic uncertainties are determined to stem from the photon energy cut, 
the impact parameter cut, fit range and signal width, description of the neutral pion 
spectrum, and the uncertainty in the f] and uj branching ratios. Systematic uncertainties 
are determined by varying cuts and taking the largest change in the calculated rate as the 
systematic error. Additionally, a systematic error is assigned representing the uncertainty 
in the extrapolation of the measured rate to Xp = 0. Details of the systematic checks are 
given below. 

A cut of 0.8 GeV is applied on the energy of photons. This cut is varied by ±0.1 
GeV. A do cut of 0.5 cm is applied to remove tracks originating from decays of particles 
at some distance away from the primary interaction point. As this cut is tight a check 



is performed to determine any associated systematic uncertainty. This cut is decreased 
to 0.3 cm and increased to 1.0 cm. Other cuts on charged track parameters do not have 
significant systematic uncertainties associated with them. 

In the fitting procedure mass spectra are fitted six times with a different fit range. 
The value for the extracted signal is seen to vary; a systematic error is assigned to this, 
calculated as the standard deviation of the six results. For the rj, significant variations in 
extracted signals occur for different choices of the background function, and a systematic 
error is assigned for this fitting uncertainty. The fitting procedure is repeated with the 
width of the signal function as a free parameter, and the variation in the number of 
extracted mesons is taken as a systematic error. To check the effect of fixing the mass in 
the fitting procedure, the procedure is repeated with the mass as a free parameter in the 
fit. The fitted masses are found to be stable and no significant variation in the extracted 
number of mesons was found, therefore no systematic error is assigned. Uncertainties in 
the meson reconstruction efficiency arise from the modelling of the 7r° spectrum. It has 
been checked that bin-to-bin migration effects are below 1% for all but the highest Xp bin. 
Even there the estimated effect is smaller than other relevant systematic uncertainties. 

For the rj, the effect of the matrix element proposed by 0] and which is not 
included in Jetset is corrected for by reweighting Jetset to agree with proposed 
parametrizations. Although little effect on the momentum spectra is seen for high t] 
momenta, the effect at low momenta is large (23% in the lowest measured momentum 
interval). As a correction, the mean of the two weighting schemes is applied to the 
measured rates, and the difference in the two schemes in each momentum interval is 
taken as the systematic error. 

Statistical and systematic errors are summarized in Tables |^ and |^ for the rj and uj 
measurements, respectively. For the systematic errors the individual errors from each 
source are shown for each measured momentum interval (1 to 6). The error for the total 
measured momentum range, calculated for each error source, is taken as the sum of the 
errors in each interval weighted by the rate in each interval. The calculations take into 
consideration whether the errors are correlated or not. The statistical errors are taken 
from the fitting procedure. The last line shows the final error for each momentum interval. 

The extrapolation of the measured rates to Xp = relies on the Monte Carlo to give 
the correct scale factor for the extrapolation. For the extrapolation the shape of the 
fragmentation function in Jetset 7.4 is used. To estimate the uncertainty in the scale 
factor the calculation is repeated for the Herwig Monte Carlo. For the t], extrapolating 
from Xp = 0.10, the difference, 5.8%, between Herwig 5.9 and Jetset 7.4 is taken as 
the systematic uncertainty in the extrapolation. For the u, extrapolating from Xp = 0.05, 
the difference, 2.4%, between Herwig 5.9 and Jetset 7.4 is taken as the systematic 
uncertainty in the extrapolation. 



5.4 Results 

Table ^ shows the results for rates and differential cross sections in each measured 
momentum interval. Results of summing over the measured Xp intervals are shown; the 
final line gives the result of extrapolating this to Xp = together with an extra error of 
5.8% and 2.4% representing the uncertainty in the extrapolation of the measured r] and 
UJ rates, respectively. Measured differential cross sections are compared to Monte Carlo 
predictions in Figs. ^ and || for the r] and u, respectively. Table ^ compares the total 
rate in the data to the Monte Carlo predictions. The errors are the quadratic sum of the 
statistical and systematic contributions. The final column gives the result of extrapolating 



Table 2: Systematic and statistical errors on the r/ production rate in each measured momentum 
interval. All values are expressed in percent. 



Source of error [rj) 


measured Xp interval 


all 


1 


2 


3 


4 


5 


6 


Photon energy cut 


4.5 




12.1 


6.7 


3.4 


1.2 


1.4 


Track selection 


1.1 


- 


2.5 


3.2 


0.4 


1.7 


1.7 


Fit range and background 


2.0 




5.0 


3.3 


3.0 


1.4 


2.0 


Fit width 


1.5 




3.2 


4.5 


0.4 


2.6 


1.4 


7r° spectrum 


2.8 




2.8 


2.8 


2.5 


2.9 


4.8 


7] branching ratio 


2.2 




2.2 


2.2 


2.2 


2.2 


2.2 


7] matrix element 


2.3 




3.9 


2.2 


1.4 


1.1 


0.3 


Total systematic error 


6.8 




14.7 


10.2 


5.8 


5.2 


6.2 


Statistical error 


3.1 




7.9 


5.6 


3.5 


2.3 


3.0 


Total error 


7.5 




16.7 


11.6 


6.8 


5.7 


6.9 



Table 3: Systematic and statistical errors on the uj production rate in each measured 
momentum interval. All values are expressed in percent. 



Source of error {uj) 


measured Xp interval 


all 


1 


2 


3 


4 


5 


6 


Photon energy cut 


4.1 


10.6 


2.2 


2.8 


1.7 


0.1 


1.8 


Track selection 


2.6 


6.7 


1.9 


1.2 


1.1 


0.5 


1.1 


Fit range 


1.0 


1.7 


3.2 


2.2 


2.0 


0.3 


1.1 


Fit width 


0.8 


1.0 


3.3 


1.8 


0.4 


0.2 


0.5 


tt" spectrum 


2.7 


2.5 


2.8 


2.8 


2.5 


2.9 


4.8 


UJ branching ratio 


0.7 


0.7 


0.7 


0.7 


0.7 


0.7 


0.7 


Total systematic error 


5.7 


13.0 


6.2 


5.1 


3.9 


3.0 


5.4 


Statistical error 


3.2 


7.8 


4.2 


3.3 


2.8 


2.0 


2.9 


Total error 


6.5 


15.2 


7.5 


6.1 


4.8 


3.6 


6.1 



the measured rate to Xp = 0. Values inside the brackets indicate the discrepancy between 
Monte Carlo and data in terms of the total experimental error. Table |^ compares the 
results of the t] measurement with results published by ALEPH [Q], OPAL and L3 |p. 
For the u, Table ^ compares the results with those of published results by ALEPH 0|, 
OPAL § and L3 §. 

For the rj, the total estimated production rate of 1.200 ±0.037stat±0.082sys±0.070extrap 
per event is somewhat higher than the predictions of both Jetset 7.4 (1.00 r] per 
event) and Herwig 5.9 (1.04 rj per event). This result is also somewhat higher than 
the results published by LEP experiments. The two measurements for the rj '-fj 
mode and the r] —>■ tt^tt^tt^ mode differ by more than twice the estimated error, the 
ratio (rate obtained with the n~^7r~7T^ mode to the rate obtained with the 77 mode) 
being 1.26 ± 0.04stat ± O.lOgys- This is consistent with a recent finding by OPAL of 
1.14 ± 0.07stat ± 0.13sys for the same ratio. 

For the u, the total estimated production rate of 0.996±0.032stat±0.057sys±0.024extrap 
per event lies significantly below the prediction of Jetset 7.4 (1.31 u per event) and 



Table 4: Measured production rates per hadronic event and differential cross sections for the 
r] and uj. The errors correspond to statistical and systematic uncertainties. The value for the 
total systematic error is taken from Table ^ and Table ^ for the r] and uj respectively, where 
correlations between errors are taken into account. The results of summing over the measured 
Xp intervals are given, including the extrapolation to the full Xp range with an additional error 
representing the uncertainty in the extrapolation. 



Xp 


1] rate 


l/o-tot ■ da/dxp 


0.10-0.15 


0.1236 ± 0.0098 ± 0.0182 


2.471 ± 0.195 ± 0.363 


0.15-0.20 


0.0753 ± 0.0042 ± 0.0077 


1.506 ± 0.084 ± 0.154 


0.20-0.30 


0.0802 ± 0.0028 ± 0.0047 


0.802 ± 0.028 ± 0.047 


0.30-0.50 


0.0611 ± 0.0014 ± 0.0032 


0.306 ± 0.007 ± 0.016 


0.50-1.00 


0.0146 ± 0.0004 ± 0.0009 


0.029 ± 0.001 ± 0.002 


0.10-1.00 


0.355 ± 0.011 ± 0.024 


all 


1.200 ± 0.037 ± 0.082 ± 0.070 




Xp 


UJ rate 


l/fTtot ■ da/dxp 


0.05-0.10 


0.2222 ± 0.0173 ± 0.0289 


AAU ± 0.347 ± 0.578 


0.10-0.15 


0.1246 ± 0.0052 ± 0.0065 


2.493 ± 0.105 ± 0.130 


0.15-0.20 


0.0825 ± 0.0027 ± 0.0042 


1.650 ± 0.054 ± 0.084 


0.20-0.30 


0.0852 ± 0.0024 ± 0.0033 


0.852 ± 0.024 ± 0.033 


0.30-0.50 


0.0569 ± 0.0011 ± 0.0017 


0.284 ± 0.006 ± 0.009 


0.50-1.00 


0.0139 ± 0.0004 ± 0.0008 


0.028 ± 0.001 ± 0.002 


0.05-1.00 


0.585 ± 0.019 ± 0.033 


all 


0.996 ± 0.032 ± 0.057 ± 0.024 
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Figure 3: Measured differential cross sections for the rj in comparison with Monte Carlo 
predictions. The errors shown are the quadratic sum of statistical and systematic contributions. 
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Figure 4: Measured differential cross sections for the uj in comparison with Monte Carlo 
predictions. The errors shown are the quadratic sum of statistical and systematic contributions. 



Table 5: Comparison of the total measured production rates with Monte Carlo predictions. 
Values in parentheses indicate the discrepancy between Monte Carlo and data in terms of the 
total experimental error. 





(r/) Xp > 0.10 


(r/) all Xp 


data 


0.355 ± 0.026 


1.20 ± 0.11 


Jetset 7.4 


0.296 (-2.3a) 


1.00 (-1.8a) 


Herwig 5.9 


0.289 (-2.5a) 


1.04 (-1.5a) 






(u;) Xp > 0.05 


(u;) all Xp 


data 


0.585 ± 0.038 


0.996 ± 0.070 


Jetset 7.4 


0.771 (+4.9a) 


1.310 (+4.5a) 


Herwig 5.9 


0.678 (+2.4a) 


1.125 (+1.8a) 



Table 6: r/ rate measured in this analysis compared to values published by LEP experiments. 



Experiment 


rj decay mode 


f] rate 


L3 i 


77 


0.93 ±0.01 ±0.09 ±0.06 


OPAL i 


77, 7r"'"7r~7r° 


0.97 ±0.03 ±0.10 ±0.04 


this study 




1.20 ±0.04 ±0.08 ±0.07 


ALEPH [| Xp > 0.10 


11 


0.282 ±0.006 ±0.015 


this study Xp > 0.10 




0.355 ±0.011 ±0.024 



Table 7: Results for the uj production rate measured by three LEP experiments. 



Experiment 


u rate 


ALEPH [0 


1.07 ±0.06 ±0.12 ±0.04 


L3 i 


1.17± 0.09 ±0.15 


OPAL |[ 


1.04 ±0.04 ±0.13 ±0.03 


this study 


1.00 ±0.03 ±0.06 ±0.02 
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Figure 5: Invariant mass distribution m{p,~^fi~) for hadronic events containing a muon pair 
after standard cuts as described in the text. The u signal of Jetset is normalized to the signal 
in the data shown. 
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Figure 6: Invariant mass distribution m(fi~^fi~) for hadronic events containing a muon pair 
with tighter cuts (as explained in the text). The u signal of Jetset is normalized to the signal 
in the data shown. 



somewhat below that of Herwig 5.9 (1.13 u per event). Published results from L3 and 
OPAL for the production rate of the u are in good agreement with the result of this 
analysis. 



6 The Branching Ratio BR(ct; — > 1^^ ) 
6.1 Extraction of the Signal 

The signal uj fi~^fi~ is obtained from the invariant mass distribution of two identified 
muons of opposite charge originating from a common vertex (Fig. ^). The signal is taken as 
the number of events in the mass region 770MeV/c^ < m(/i"'"/i~) < 795MeV/c^ reduced 
by the estimated background in the same region. Monte Carlo studies have shown the 
background to be linear as a function of the invariant mass. The level of background in 
the signal region is extrapolated from the background in the sidebands in the mass ranges 
500MeV/c2 < m(/i+/i-) < 745MeV/c2 and 820MeV/c2 < m(/i+/i-) < lOOOMeV/c^. 
Applying the cuts described in Sections ^ and ^ results in an optimal signal-to-background 
ratio, S/ y/S + B. 



6.2 Measurement of BR(a; — ^ M"*"/^ ) 

Figure ^ shows 35 events in the signal region. A background of 16.94 events is estimated 
from the sidebands, leaving 18.06 ± 5.92 signal events. When the acceptance is taken 
from Jetset 7.4, but with the momentum spectrum and rate corrected to agree with 
the 7r+7r~7r'^ analysis, a muonic branching ratio BR(ci; fJ'^fJ' ) = (9.0 ± 2.9) x 10~^ is 
obtained, based on the selected events shown in Fig. ^ 

The effect of tightening cuts in order to increase S/ \/B is shown in Fig. ^ For this, 
the additional requirements are muon momenta greater than 2.5GeV/c, two VDET hits 
instead of one, and a stronger bb event rejection (keeping 90% of light quark-antiquark 
events). In the signal region 14 events are kept, with an estimated background of 3.24 
events, leaving 10.76 ± 3.74 signal events. 

6.3 Systematic Error Analysis 

Uncertainties on the muonic branching ratio come from the uncertainty of the momentum 
spectrum of the u and its rate in hadronic events, which are taken from the tt+tt^tt^ 
analysis. If the momentum spectrum from Jetset 7.4 is used, the extracted branching 
ratio would increase. Half of the difference from the shape of the spectrum is taken 
as systematic error (4%). Further systematic errors come from the uncertainty of the 
measured rate of a; production per event in the u — > 7r+7r~7r° analysis. Because event 
selection and track selection are similar in both analyses their errors are assumed fully 
correlated and neglected. The remaining errors (vr^ selection, fitting procedure, etc.) are 
studied by varying the Monte Carlo momentum spectrum according to these remaining 
errors. An error of 6% is obtained. 

Possible uncertainties can arise from the track selection and b-tag probability. The 
last also changes the flavour composition of the sample; the u rate for different flavours 
has not yet been measured. The tt+tt^tt'^ analysis was therefore repeated with the cuts 
on b-tag probability and track selection from this analysis. An error of 4% is derived. 
The estimation of the background is dominated by the limited statistics in the sidebands. 
Smaller uncertainties come from a possible curvature of the background shape, e.g., due 
to the p meson. The error is determined to be 6%. Further systematic uncertainties are 
related to the performance of the ALEPH detector [|1^: muon identification (±5%), 
tracking resolution, which influences the signal shape (±5%), and the missing energy, the 
impact of which is found to be negligible. The limited Monte Carlo statistics adds 3% 
to the error. Other sources of systematic error such as background are found to be small 
and are neglected. The total systematic error is obtained by adding the contributions in 
quadrature, giving ±13%. 

6.4 Results 

The muonic branching ratio is measured to be BR((X' f^^f^ ) = (9-0 ± 2.9 ±1.1) x 
10~^ from the observation of 18.06 ± 5.92 signal events. This is in agreement with 
the expectation from the flavour composition of the vector mesons and from lepton 
universality. 



7 Summary and Conclusion 



The inclusive production of the t] meson and the ui vector meson in hadronic Z decays has 
been studied and compared to model predictions. Decays of the type r] tt+tt^tt^ are 
reconstructed for the momentum interval of Xp > 0.10 where Xp = Pmeson/pbeam- Decays 
of the type u — > 7r"'"7r~7r° are reconstructed for the momentum interval of Xp > 0.05. A 
signal is seen in the muonic decay u — > ■ 

The average rj rate per event for Xp > 0.10 is measured to be 0.355±0.011stat±0.024sys; 
an extrapolation to Xp = using the shape of the fragmentation function in Jetset 7.4 
yields a total production rate of 1.200 ± 0.037stat ± 0.082sys ± 0.070extrap per event. The 
predictions of both Jetset 7.4 (1.00 rj per event) and Herwig 5.9 (1.04 rj per event) 
are somewhat below this result. The use of the matrix elements proposed by ^ and 
20| resulted in a significant effect on the momentum spectrum for low rj momenta. The 



matrix element is not present in Jetset 7.4 and Herwig 5.9. Published results of L3 
0] and OPAL |^ are somewhat lower than the result of this analysis. 

The average u rate per event for Xp > 0.05 is measured to be 0.585±0.019stat±0.033sys; 
an extrapolation to Xp = yields a total production rate of 0.996 it 0.032stat ± 0.057sys ± 
0.024extrap pcr eveut. The rate lies significantly below the prediction of Jetset 7.4 (1.31 
uj per event) and somewhat below that of Herwig 5.9 (1.13 u per event). Published 
results from L3 and OPAL for the production rate of the u are in good agreement with 
the result of this analysis. This measurement also agrees with that previously published 
in PI and improves it substantially. 

In the muonic decay mode of the u vector meson, 18.1 ± 5.9 events are observed 
yielding the first measurement of the muonic branching ratio of BR(ti; —>■ /i"*"/!") = 
(9.0 ± 2.9stat ± 1-lsys) X 10~^. The result is in agreement with the expectation from 
the flavour composition of the vector mesons and from lepton universality. 
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